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ON-BOARD DATA MANAGEMENT STUDY FOR EOPAP 


1. Summary 

A study of the requirements, implementation techniques and mission 
analysis associated with on-board data management for EOPAP has been 
performed. SEASAT-A has been used as a baseline with reference M 
used as the primary source of information for the spacecraft config- 
uration. The storage requirements, data rates, and information extrac- 
tion requirements have been investigated for each of the following 
proposed SEASAT sensors: a short pulse 13.9 GHz radar, a long pulse 
13.9 GHz radar, a synthetic aperture radar, a multispectral passive 
microwave radiometer facility and an infrared/visible very high resol- 
ution radiometer (VHRR) . The short pulse radar is used for altimetry and 
the determination of wave height and wave directional spectra. Primary 
consideration in this study was given to this sensor based on data pro- 
vided by the GSFC technical officer. The long pulse radar is used for 
measuring wind speeds by scatterometry . The synthetic aperture radar 
provides side-looking coherent imagery. The passive microwave sensors are 
radiometers at 5,13.9,18,22,36, and 53 GHz which are used for correcting 
altimeter measurements, measuring surface roughness and sea ice. The 
VHRR sensor provides cloud cover maps and sea surface temperature infor- 
mation. Data management principles for the latter were included in the 
study based upon the results of the earlier study contract NAS 5-21940 
"ITOS VHRR On-Board Data Compression Study" (reference ^2j) and extended 
to include the other SEASAT sensors. 

Rate distortion theory was applied to determine theoretical minimum 
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data rates and compared with the. rates required by practical techniques. 
It was concluded that practical techniques can be used which approach the 
theoretically optimum based upon an empirically determined source random 
process model. 

Finally the results of the preceding investigations were used to 
recommend an on-board data management system. The purpose of the on-board 
data management system is for (1) data compression through information 
extraction, optimal noiseless coding, source coding with distortion, data 
buffering and data selection under command or as a function of data act- 
ivity, (2) for command handling, (3) for spacecraft operation and control, 
and (4) for experiment operation and monitoring, it was concluded that 
on-board general purpose programmable computer should be used and that 
the advanced on-board processor (AOP) developed by NASA GSFC is suitable 
for the tasks required. 

Section 2 is a review of the SEASAT sensors studied and their data 
rate requirements. Section 3 presents rate distortion theoretic consid- 
erations. Section 4 is a discussion of the data reduction and estimation 
techniques applied to aircraft samples of the short pulse radar type and 
the computer programs developed in this study for the analysis. Section 
5 presents results for the VHRR data. Finally section 6 discusses an on- 
board data management system. 


\ 
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2. . SEASAT Sensor Data Processing 
2.1 Short Pulse Radar 

The short pulse radar is used to measure the mean altitude of the 
spacecraft and the rms wave height in a nadir-looking mode with a .6m 
dish. It also can be used with a 2 m dish in an off-nadir conical scan 
mode to measure wave directional spectra. The radar uses linear chirp 
FM transmission to attain the required resolution («3 nsec compressed 
pulse length) in altitude, wave height, and wave directional spectra. On- 

J 

board processing can be used to estimate and transmit these parameters 
at a low bit rate relative to the bit rate required for the unprocessed 
received radar data. The actual bit rate resulting depends on the number 
of read-outs per second of the processed estimates which can be control- 
led by the on-board data management system. The control can be through 
read-out rate selection under ground station control or under data con- 
trol as a function of data activity. In either case, data compression 
can', be used for further bit rate reduction. by sending encoded changes 
in the parameters rather than the parameters themselves. In this method 
noseless coding . is used to take advantage of the lowered entropy of the 
difference sequence. 

Further details regarding the processing of the radar data for 
parameter estimation appears in section 4 based upon samples of aircraft 
data supplied by the NASA GSFC contract technical officer. 
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2.2 Long Pulse Radar 

The long pulse 13.9 GHz radar is used as a scatterometer to measure 
wind speed and direction. The scatterometer shares the 2 m dish with the 
short pulse radar in a conical scan mode and operates on the principle 
that the radar cross section of the sea surface depends on the wind 
speed and direction. Atmospheric calibration is provided by the 13.9 GHz 
passive radiometer. As for the short pulse radar, on-board processing can 
be used to reduce the received radar data into estimates of wind speed 
and direction so that arbitrarily low bit transmission rates can be 
attained depending upon the read out rate desired, and/or the data ac- 
tivity and/or the use of data compression on the sequence of parameter 
estimates. 

2.3 Synthetic Aperture Radar 

A synthetic aperture radar is provided for side- looking coherent 
imaging. Waves on the ocean can be analyzed from the imagery to determine 
wave spacing and direction (except for a 180° ambiguity) . From the wave 
periodicity the rms wave height can be found through the known determin- 
istic relation between the two parameters. Furthermore information on. 
icebergs, shoals, kelp beds, tides, oil slicks, etc. can be provided. The 
major drawback in the useage of this technique is the high data rate. 

This data rate can be reduced by selective coverage through the on-board 
processor and through prefilitering and buffering operations to 8 Mbps. 
Because of the high data rate and the inherent redundancy of images, 
particularly the synthetic aperture type , the signal is an ideal candi- 
date for data compression. Using rate distortion theory (see section 3) , 



5 


and empirical experience with images, it is estimated that the bit rate 
could be reduced to 1 Mbps with acceptable quality by using a combina- 
tion of line-to-line and sample-to- sample compression techniques. If 
further degradation is allowed so that only the essential information is 
retained (wave locations, edges of land or iceberg or other significant 
formations, etc) a reduction to .5 Mbps or less may be possible. 

2.4 Passive Microwave Radiometer Facility 

The passive microwave radiometers provide atmospheric path length 
corrections for the short pulse radar altimeter, sea surface temperature 
measurements, surface wind speed measurements, and rain area and cloud 
water content determinations. The bit rate requirement for this data is 
3.7 kbps. This can be substantially reduced by data compression using 
the inherent redundancy of the data (see section 3) . 

2.5 Very High Resolution Radiometer 

A two channel very high resolution infrared/visible radiometer (VHRR) 
provides cloud cover pictures and surface temperature measurements in the 

10.5 to 12.5 micrometer infrared band and in the 0.6 to 0.7 micrometer 
visible band. This instrument is essentially identical to that used for 
ITOS-D. The scan voltage is a 35 kHz low pass signal which requires on 
the order of 1 Mbps for full-resolution, uncoded transmission. However, 
under samp ling, line deletion, data formatting, and data compression 
under control of the on-board processor can be used to substantially 
reduce the required bit rate. Detailed methods of compression for this 
source data appear in Systems Analysis report no. 75100, "ITOS VHRR On- 
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Board Data Compression” prepared under NASA GSFC contract NAS5-21940. A 
summary of the essential results appears in section 5 of this report. 
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3. Rate Distortion Theory 

The rate-distortion function of a source with a given probability 
distribution determines the minimum channel capacity required to transmit 
the source output (or the storage capacity required to store the source 
output per unit time) as a function of the desired maximum average 
distortion, where the distortion function (sometimes called fidelity 
criterion) is a measure of the agreement between the source output 
and the end user output as specified by the the user. The concept ori- 
ginated with Shannon who calculated the rate-distortion function for 
certain sources including in particular a Gaussian discrete time random 
process. This was later extended to more complicated sources and applied 
to performance bounding by other authors. 

Figure 3.1 presents a highly simplified communication system block 
diagram. The system designer is given a source and wishes to encode the 
source in such a way that the channel capacity (or storage) require- 
ment is minimized. To make this minimum as small as possible, he is 
willing to tolerate some average distortion (typically average error pow- 
er) between the source output and the decoder output. The problem ad- 
dressed by rate distortion theory is the minimization of the channel 
capacity requirement while holding the average distortion at or below an 
acceptable level. 

To be slightly more specific, in figure 3.1 let the average infor- 
mation transmitted from the source to the decoder output be denoted by 
the function I(X,Y), more information corresponding to larger values of 
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I(X,Y). Let the analogous average information transmitted from the enco- 
der output to the channel output be denoted I(U,V). The channel capacity 
C is defined as the maximum of I(U,V) over all possible input devices. It 
is obvious that in figure 3.1, the intervening nature of U and V implies 
I (X,Y) <1 (U,V) <C if the function !(*,•) is properly defined. Shannon 
defined the correct function for information and defined the rate 
distortion function R(D) as the minimum value of I(X,Y) for a given 
distortion level' D. He also showed that R(D)<C issures the possibility 
of attaining average distortion D by coding at rate C with a sufficiently 
complicated encoder . 

To look at the problem from a different angle, the encoding 
and decoding can be separated into two parts -source coding and channel 
coding. The job of the source coder is to remove source redundancy, where- 
as the job of the channel coder is to insert controlled redundancy to 
combat noise. How well channel coding can perform is established by 
Shannon's noisy coding theorem which states that transmission without 
channel errors is possible if and only if the rate at which messages are 
presented to the channel coder from the source coder does not exceed 
channel capacity. With this constraint on the source coder, at a given 
channel capacity the source coding must result in a loss of quality-i.e. 
some distortion must result. The rate distortion function R(D) is the 
minimum source coder rate (minimum required channel capacity) for an 
average distortion level D. 

An analytically convenient source model to use is the Gaussian which 
frequently is a good model in practice, in particular for the image 
data found on SEASAT. Consider the problem of encoding a source which 
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generates a sequence of statistically independent Gaussian random var- 
iables with a squared error distortion measure- This source is char- 
acterized by a known mean (i.e. dc) value, which does not affect coding 

2 . 

rate because it is a simple shift, and a variance (power) a . Let the 

source output be the N values {x^,i=l,2 , - - . -N} and the decoder output 

be (y^,i=l,2, .N>. The distortion is taken as 

N 2 
d(X,Y)= Z (x i -y i ) 

i=l 


so that the average distortion is the average mean square error (error 
power) : 

1 N 2 
D =i SEUi-yi) 

.. i-1 

It can be shown that the rate distortion function is given by 

2 


2 /? 


R<D) = 


ijlog^- = a > 


a 2 £ D 


(3.1) 


That is, the rate distortion function is one-half the logarithm of the 
signal power-to-distortion (noise) ratio. It is achieved conceptually 
by encoding in such a way that the output error is Gaussian with variance 
D on each sample value and is sample -to -sample independent, but not 
necessarily .signal independent. Actually there are no known practical 
methods' for achieving a Gaussian error distribution or for achieving 
the rate distortion bound of (3.1). Instead one must be content with 
quantization methods which can be made to be within approximately 
.25 of a bit or so of the bound by choosing the quantization levels 
properly and using variable length coding on the quantized values. 
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In terms of quantizing schemes the rate distortion function has 
the following rough intuitive justification. The noise standard devia- 
tion as a fraction of the signal amplitude is inversely proportional to 
the number of quantization levels. Therefore, the number of levels should 
be proportional to cr//B and the number of information bits should there- 
fore be the logarithm (base 2) of this quantity. This is what is stated 
in (3.1). If the allowed distortion is greater than the signal power, i.e. 

D >0 , the minimum transmission rate and hence rate distortion function 
must be zero since nothing need be transmitted at all to achieve error 
power equal to signal power. This is also stated in (3.1). 

The preceding was for a source of independent values. Suppose now 
the source contains memory as in most physical sources. In particular, 
suppose the source output is Markov (which means roughly that the future 
depends only on present values and not the past) with correlation coef- 
ficient p , jp| <1. This is a model which has found good agreement in practice, 
in particular for video data of the VHRR and synthetic aperture radar type. 
The value p is the correlation between adjacent values on a scan line. 

If p= 1, then each value is the same as the last. If p =0, the values 
are independent as before. For this model it can be shown that (3.1) 

for. low distortion levels becomes 

2 2 

R(D) = h log ° ‘ (3.2) 

The typical range of the parameter p is .95 to .99. Thus the 

2 

number of bits can be reduced by ^ log(l-p ) , or something up to 3 bits, 
for a given mean square error D over coding for independent values. This 
perfprmance can be approached through differential coding methods whereby 
successive sample differences are encoded by a variable length code. For 
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image data further gains can be made by using line-to-line effects. 
Assuming the same correlation model between lines, another gain of up to 
3 bits can be attained assuming sufficient random access storage is 
available to store one whole line of data . (in addition to the 
present line ) . 

Rate distortion theory is most useful in practice in providing 
a bound on the performance of practical data compression methods. This 
is the primary way it has been used in this study. Consider for example 
differential variable length coding wherein the successive sample 
differences are quantized and encoded by a variable length code. Assuming 
fine quantization, the variance of the difference between the present 
quantized value and the next sample is approximately 

2 (l-p);»<l+p) (1-p) ** 1-P 2 (3.3) 

Allowing .25 bit for the optimum quantized value code over the best 
that can be done as shown by (3.1) and (3.2), it is seen that the dif- 
ferential method comes within the .25 bit of the optimum up to the approx- 
imation of fine quantization and the approximations of (3.3) 
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4: Short Pulse Radar Data Reduction 


4.1 Introduction 

Extensive studies were performed on the extraction of rms wave 
height and wave periodicity based on off-nadir looking aircraft radar 
measurements provided by the NASA GSFC contract technical officer on 
digital 9 track tape. Two aircraft runs were provided, the first being 
from the JONSWOP experiment with the NRL altimeter used in an off-nadir 
mode by banking the aircraft. In this experiment a 10 nsec pulse width 
was used with a 6° beamwidth antenna. The aircraft was banked so that 
the beam center pointed about 10° and 15° off-nadir on passes in various 
turning configurations. The aircraft was at an altitude of approximately 
11,000 feet. Similar data was provided by an NRL experiment in 

January 1974, the primary differences being a doubled sampling rate, a 
fixed bearing . rather than a turning course by the aircraft, and less 
noisy data. 

The sampled data provided in each of the experiments consisted 
of 6 digital tape files corresponding to each of 6 aircraft passes under 
different conditions. Each pass consisted of 8192 signal traces measur- 
ed by sampling the pulse returns sequentially 160 times with increasing 
delay on each pulse across a total sampled time interval of one micro- 
second for the JONSWOP data and 500 nsec for the NRL data. A typical pulse 
appears in figure 4.1.1 and is seen to comprise the total returned 
energy from the sea surface plus noise with noise only preceding and 
following the return signal. 
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Computer programs were developed for the analysis of the sampled 
data and are included in listed form as a portion of this report. These 
are in Interdata assembly language and FORTRAN IV. The analysis programs 
include an Interdata assembly language program for inputting the data 
from the digital tape under a FORTRAN call and a master program with the 
ability under input control to print out individual signal traces or 
averages of signal traces, the ability to print out the Fourier transform 
amplitude spectra of the signal traces or their averages or the average 
of the Fourier amplitude spectra of the individual signal traces, using a 
fast Fourier transform routine writen for the Interdata in FORTRAN, the 
ability to print out correlation functions and cross-correlation functions 
of the signal traces, and the ability to print out histograms of the 
signal trace peaks subject to a threshold on the signal level, includ- 
ing the calculated mean and rms values of the histogram. Except for 
the assembly language digital data input routine, FORTRAN IV was 

used exclusively so that the programs could be compiled on any machine, 
except, possibly for some minor modifications to conform with the 
particular compiler conventions . 

Extensive computer program outputs of the various types were pro- 
vided to NASA GSFC during the course of the contract for the purposes of 
analyzing the various approaches to determining the wave heights and/or 
wave periodicities from off— nadir short pulse radar data. The results 
to date appear promising but inconclusive, partly due to the limited 
amount of data, the noisiness of the data, the lack of resolution of the 
data, the necessity to construct one pulse from 160 sampled pulses, and 
the lack of ground truth values. 
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4.2 Simplified Theory of Wave Height and Wave. Periodicity Estimation 
M^ny papers have been written on the theory of wave height and 
wave period estimation from radar sea scatter measurements. Some of 
these are listed in the references. In this section a brief review of 
of the pertinent theory to off— nadir short pulse radar measurements is 
given in highly simplified form. 

Consider a radar illuminating a periodic sea surface at an angle 
a from nadir as shown below: 

radar 



Figure 4.2.1. Short Pulse Radar Geometry 

The returned signal power is the radar cross section as a function of 
propagation path length across the illuminated surface. If one assumes 
that the returns move linearly across the surface in time, the returned 
signal power contains an amplitude modulation component of frequency 
c/(2Lsina). This conceptually can be extracted either directly from the 
returned power as a function of time or from the Fourier transform (spec- 
trum) of the returned power by peak determination. Complicating the task, 
of course, is the non-pure periodicity of the sea surface, the nonlin- 
earity of the range-time scale across the sea surface and the receiver 
noise. In the time domain the estimate could be extracted by a bank of 
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matched filters, incorporating the range-time base correction as a 
function of the off -nadir angle, a . Alternately the data could be samp- 
led non-uniformly or interpolated to correct the time base followed by a 
Fourier transform of the squared envelope. The peak of the spectrum 
(away from D.C.) represents the estimated amplitude modulation frequen- 
cy. 

Given the wave period estimate, the rms wave height estimate can 
be determined from the deterministic relation between the two. Alternately 
the wave height can be estimated and used to estimate the wave period 
through the same relationship. The average power return over a number of 
pulses is just the convolution of the transmitted pulse with the prob 
ability density function of the wave height integrated over range. The 
spectrum is the product of the transform of these factors. Thus the rms 
wave height can be extracted from the returned signal by finding the rms 
value of the probability density imbedded in the averaged returns. Ref- 
erence [l\ proposes doing this from the leading edge of the return which 
should be more or less spread out depending on wave height (actually [l] 
proposes using nadir measurements although the same analysis can be ap- 
plied to small off-nadir observations). Levine [3,4] shows that the spec- 
trum should be approximately proportional to 

t 2iro 2 1 2 tt £ 

2 ( - - cosav ) J sine (— ^— sinct ) 

where : 

V = frequency 

I(v)— Fourier transform of the transmitted pulse squared 

a — rms wave height 
z 
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a= off-nadir angle 

illuminated surface length. 

This is just the product of the transforms of the above referred to 
factors assuming the wave height distribution is Gaussian. Thus can 
be estimated by the rate of exponential decay of the spectrum. 

Alternately, Levine [V] has shown that 0^ can be estimated from 
a histogram of the locations of peaks in the returned signal. In par- 
ticular, he found that 

2 1. 2 U U sina) 2 

a z = sec a jjc/2) r J 

2 

where T is the rms peak arrival time. Thus a could be extracted from 
n z 

a peak histogram over several pulses. 

Both of the preceding methods of extracting are theoretical 
results from a 2 dimensional model and involve approximations. Thus 
experimental results are needed to confirm the theory or point to 


needed modifications. 
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4.3 Computer Programs Developed for the Analysis 

4.3.1 Data Input Program 

The data as supplied by the NASA GSFC technical officer is for- 
matted on digital tape in 6 files of 257 records. The first record is a 
short header record for the file. Each succeeding record consists of 
32 signal traces, formatted with 14 header bytes followed by 160 2-byte 
signal sample values. Thus one record consists of 32 x 167 = 5344 
2-byte (halfword) values. Each file of 256 data records comprises 
one aircraft pass. 

The program written to read these values into a FORTRAN array 
is writen in Interdata assembly language and appears as appendix A. The 
program is called by: 

CALL ECKERM(INP) 

where INP (N) is a (2-byte) integer array of dimension not less than 5344. 
All 5344 values from the next record on tape are unpacked sequentially in 
the array INP, including the header values. 

4.3.2 Fast Fourier Transform 

A fast Fourier transform was used extensively for data analysis 
and is listed in Appendix B. It was written in Interdata FORTRAN IV and 
is a subroutine called by: * 

CALL FOURlT (DATA,N,S ,ISIGN) 

where : 

DATA (I) = complex array of N values 

S = table of sines of dimension not less than N/2 + 1 



ISIGN is used as follows: 
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sign (ISIGN) = sign of the exponent in the Fourier transform 
If [isign| f 1, the sine tables are set up in S 
If ISIGN = 0, tables are set up only. 

The quantity calculated is : 

N 

E DATA(k) exp (i -4-(k^.l) (n-1); sign (ISIGN)) , n=l,2,---,N 

k=l 

The results appear in DATA (n) in place of the original data. Note that 
no normalization by N occurs. 

4.3.3 Analysis Program 

The listing of the main analysis program appears in Appendix C as 
it is currently configured. Many modif ications were made to the program 
as the contract progressed to conform to the various analytical approaches 
and data requests by NASA GSFC. Certain classes of outputs can be gen- 
erated by an appropriate input to the program from the teletype 
keyboard as the program now stands. Others can be generated by simple 
program modifications and re-compilitation. The following outputs have 
been generated to date: 

(1) . A printout of individual signal traces squared (i.e. power) 

or the average of a number of signal traces squared, with 

« 

the number under input control , and the centroid of the 
printed quantity. An example of this printout appears in 
Appendix D. 

(2) . A printout of the log of the spectrum of the average in (1) . 

An example appears in Appendix D where a 1024 point spectrum 
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appears, generated by setting all values to zero outside 
the return signal range. Fourier transform amplitude spectra 
can also be generated without the appended zeroes, in which case 
the values are windowed prior to transforming. 

(3) . The average of the Fourier transform amplitude spectra of 

individual pulses . 

(4) . Correlation functions of the signal traces and cross corr- 

elation functions between signal traces. 

(5) . Numbers (2) -(4) with the range-time base corrected by inter- 

„ polation to be linear. 

(6) . Signal-to-noise ratio estimates. 

(7) . Histograms of peaks exceeding an input threshold, including 

the mean and rms values of the histograms. An example of 
this printout appears in Appendix E. 

To date the efficacy of the various estimation techniques remains 
unproven, although the results are promising. Upon the attainment of 
further data, the programs developed under this contract will be available 
for the analysis and the developement of an optimal technique. 
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5. VHRR Data Compression 

The VHRR signal appears as in figure 5.1 with the timing as 
indicated in table 5.1. The voltage range is 6 volts which is contained 
within a lowpass bandwidth of 35 kHz nominal. Using a sampling rate on 
the order of 100 kHz (1.43 times the Nyquist rate) and something like 8 
bit sample quantization results in a bit rate on the order of 1 Mhz for 
a straightforward PCM transmission system. However, an examination of 
table 5.1 reveals that 40 % of the data can be essentially eliminated by 
transmitting the temperature and calibration values as single digital 
values rather than as long repetitive steps and by synchronizing on and 
selecting out the earth scan data only. It is shown in the results of an 
earlier study [2^ that another 3 bits can be removed from the 8 bit 
quantized values in the earth scan data by using variable length coding 
on the sample differences as described in section 3 of this report. Thus 
the compressed-buffered bit stream requires a channel capacity on the 
order of (8-5) bits x 60% x 100kHz & 300 kHz. Further details can be 
found in the reference [2], 

As indicated in section 3, a further reduction in bit rate can be 
achieved by using line-to-line effects, perhaps to as low as 100 kHz if 
enough storage is allocated to hold the last line. Still further 
reduction in the average rate can be attained by selecting out only the 
intervals of time of particular interest or by deleting samples and/or 
lines on a regular basis. In this case digital filtering should be used 
to prevent aliasing on reconstruction. 
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Hot clock 
controlled 

43.32t 

Amplitude varies over range OX to 
100X (OX Hot, 100X Cold). Period, 
varies from 43.68 msec @900 n.ml. 
to 48,64 msec @600 n.ml. 
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16.21 v 
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Two level signal starting at OX 
amplitude for ^ the period when » 
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remainder of the period 

2n4 Sub-Sync starker 

13056 

43.51 

; i6 ./• 

0.05 
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Space scan (post earth) 

N 
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Controlled 

*52.52 

mm 

4.64t 

Amplitude Is undefined as signal 
may be influenced by the visible 
channel calibration target 

Voltage Calibration 

17152 

57.17 ; 

. ‘ 1792 

3.97 

Seven levels increasing In ampli- 
tude from 0 to 100X. Zero level 
following 7th step. 

Visible calibration 

18944 

63.14 

» 

256 

0.85 

Normal range when target la illu- 
minated by sun is about 20X amplitude 

Radiance calibration 
300k target 

-"V . 1 1 

19200 

63.99 

256 

0.85 

Amplitude about 20X. Varies with 
target's actual temperature 


TABLE T5,i'-' VHRR TIMING 


24 





After Start of IB Precursor 

Duration 

Signal Character la tics 


function 

Counts * 

Tima (msec) 

Counts 

Tine (nsec} 

• 

Temperature 8 ana or 
300K target (expanded 
range +20 to +34°C) 

19456 

64 .83 

256 

0.85 

for normal temperature range ampli- 
tude will be between 0 end 80%. 

(If the temperature la out of this 
range , the signal will either be 
0 or. between 85% and 100%) 


Temperature Sensor 

19712 

65.70 

256 

0.83 

Same as above 


(expanded) range +7 
to +21 C) 



■ "• V . 




Temperature Sensor (ax- 

19968 

66.55 -Vi, 

256 / 

0.83 

Same as above 


pended range -6 to +8°C) 



" • ■*' '' 

t 



Post-earth Synchronisation 

’ % ■ '* " 


' , .1 

i 

^ . - i 

- . . ■ . • , : ! 


front Porch 

20224 

67.41 

1024 , 

3.41 . 

Amplitude 0% 


Sync Pulaa 

21248 

. 70.82 

236 

. '0.B5 

Amplitude 100% / . 


Back Porch 

21504 

71.67 

- 256 

0.85 

Amplitude 0% , , , 

■ ( 


Ind of IB scan 

21760 ' 

72.53 



; ' . i 

1 

Precursor Via 

22523* , 

\ 

75.00+1.00 

* 512 

1.71 

t 

12.50kHi equate wave. 6 to 100% 
amplitude. Starts 72+10 before nadir 1 




’ . a 


1 • 

1 

o S 

Pra- earth Synchronisation 



V. , 

• f 

Signal 


. 



‘ i 

o.g 

Front Porch 

2303S 

76.78 ■ 

236 

0.85 

Sams as IB signal 

* * 

' 

Sync Pulse 

23291 . 

77.63 

312 , 

1.71 

Same as 1& signal - 


Back Porch 

< 23803 

79.34 

256 

0.85 

Same as 1& signal j 

.3* 

Space scan (pra-aarth) 

24059 

80.19 

' . Mot clock 
, controlled 

2.08t 

. . * 

Amplitude 0%. Period varies 
1 identically with XK signal 

♦ 


ro 

tn 


TABLE ^5^1- (cent) 


i 


/ 



ORIGINAL PAGE IS 
OF j>OOR QUALITY 


Function 

After Start 
Counts* 

of IR Precursor 
Time (msec) 

Duration 

Counts Time (msec) 

Signal Charecterlstlce 

Earth scan 

Not clock 
controlled 

82,27+ 

* , Not clock 
. controlled 

45.32+ 

Amplitude varies over range of OX to 
100X. (OX Black. 100X White). Period 
varies Identically with IR channel. 

1st Sub-Sync marker 

27387 

91.28 

16 

0.05 

Same as IR signal 

2nd Sub-Sync marker 

; 35579 

118.50 


0.05 

Same as IR signal 

Space scan (post earth) 

Not clock 
controlled 

127.59 

Not clock 
controlled 

4.64+ 

Sams as IR signal 

Voltage calibration 

39675 

132.24 

1792 

5.97 

Same as IR signal 

Visible calib, target 

41467 

138.22 

256 

0.85 , 

Same as IR signal 

Radiance calibration 
300°K target 

41723 

139.00 

* J 

236 

0.85 

Seme as IR signal 

Temperature sensor 
(+20 to 34°C) 

41979 

139,93 

f\ V- 256 : 

• 0.85 

Seme as IR signal 

Temperature sensor 
(+7 to +21°C) 

42235 

140.78 : .y. 

, tn ” 6 ‘ ■ 

•♦***. * 

0.85 

Same as IR signal " . ■* ; , . 

Temperature sensor 
(-6 to +8°C) 

42491 

141.64 

256 

0.85 

Same as IR signal . 

Foac-earth Synchronisation 8ignal • 

1 

r #• 

4 ' 

( 


Front Porch 

42747 

142.49 

1024 

3.41 

Seme ee IR signal 

Sync Pulse 

43771 

145.90 

256 

0.85 

Setae ee IR algnel 

Beck Porch 

44027 

146.76 

256 

4 

. 0.85 

Seme ee IR algnel 

Multiplex tolerance 
safety tone 

44283 

147.61 


•-2.39 


End of Scan • 


150.00 



* ■ ' V 


* A count la ona cycle of 300-kHe iqutn wave. \ 

+ Spacecraft at 790 n.ml., 0° roll error, and nominal position for start of precursor 

t Nominal value for switchover from IR to Visible Scan. All following times assume that this switchover is nominal. 


i 


TABLE 4 . 5..1 


(cont) 
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6, On-Board Data Management System 

It is clear from the preceding that the on-board processing requires 
a great deal of flexibility and computation which can be best performed 
by an on-board programmable computer , the programs for which can be . 
loaded under ground station control. Such programs can be used to select 
and multiplex the sensors, to reduce the data to the desired parameter 
estimates, select sampling rates, encode the data and buffer it for 
transmission and/or storage at a constant rate from variable rate 
input values. 

The primary considerations in the selection of an owboard pro- 
cessor are the computation speed, instruction repertoire and storage 
capacity. Also important, of course, are the power consumption, weight 
and volume. An investigation of the NASA GSFC-developed advanced on-board 
processor (AOP) indicates that it is suitable for SEASAT. It is a 1 
microsecond cycle time machine with 16K 2 microsecond cycle time memory 
banks. There are 11 load/store and 9 arithmetic instructions plus various 
I/O, interrupt, branching, and test instructions. 
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APPENDIX A 


Data Input Program 

..'READ ECKERMAN TAPE INTO. FORTRAN 5X74 


PAGE 


0000 R 


ENTRY 

ECKERM 


0000R 


EXTRN 

* 

* CALL IS* 

. 0 




X CALL ECKERMC IMP ) 




X INPCN) IS h 

FIX PT ARRAY 




X N.GE. 5344 
X 

X LDD 5/74 





X MOD FOR 2 BYTE 10/74 


0000R 

D0A0 

00S8R 

* 

ECKE'RM STM 

<10 , REG 

SAVE USER REGS 

00Q4R 

48AF 

0000 

LH 

10,0(15) 

CK NUMB ARCS 

0008 R 

27h4 

SIS 

10,4 


. 000AR 

2337 

EZS 

OK 


000CR 

C8B0 

3333 

LH I 

1 1 , C ' 3 3 ' 

SEND ERR MESS 

0010R 

41F0 

0000F 

BAL 

1 5 , . 0 

■00 1 4 R 

E i 30 
0000 

SVC 

3,0 

& STOP 

V®1 SR 

4SAF 

OK LH 

10,2(15) 

GET A( INP) 

f : y : 

.0002 



SAVE FOR LATER 

. 00 ICR 

08DA 

LHR 

13,10 

00 TER 

4 0 A0 
009 AR 

STH 

10, INPUT+4 . 

SET ’UP SVC 

, V 00S2R 

C AA0 
29BF 

. AH I 

10,10687 

END ADDR 

. 0 .0 2 SR 

40A0 

009GR 

’ STH 

10 , INPUT+6 


002 AR 

2 4E2 

LIS 

14,2 

LAST ADDR 

;.002CR 

C8FD 
£9 BE 

LH I 

15, 10686(13) 

; 0030R 

El 10 

GO SVC 

•1 , INPUT 

' READ TAPE IN 


0096R 



■CK STAT 

" 00 3 4 R 

1. ’ * .v‘ 

4800 
0 0 9 8 R 

LH 

0, INPUT+2 


0038R 

4 230 


BNZ 

NO GO 




003CR 

0064R 

D3BD 

LOOP 

LB 

11,0(13) 

; GET MSP 



0040R 

0000 

C7B0 

* 

XHI 

1 i , X ' 1 0 ' 

■ .MAKE ALL 

+ 


00 4 4R 

0010 

C4B0 


NHI 

11 ,X'1F' 

MASK OUT 

INS 

BIT 

0 0 4 8 R 

00 IF 

D3AD 


LB 

10,1(13) 

" GET LSP 



004CR 

0001 

C4A0 


NHI 

1 0 , X ' 1 F ' 

• MASK OUT 

INS 

BIT 

0OS0R 
0052R 
0 0 5 4 R 

001F 
9 1 B5 
©6BA 
40BD 


SLLS 

OHR 

STH 

11,5 

11,10 

11,0(13) 

: iLINE UP 

■ PACK. 

STORE ALJAY 


0058R 

0000 
C 1 D 0 


BXLE 

13 , LOOP 

GET NEXT 
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READ ECKERMAN TAPE INTO 1 ' FORTRAN S/?4 


003CR 


005OR 

D1A0 


LM 

1 0 , REG 

RETURN IF DONE 


0088R 





0 0 <b 0 R 

4 30F 
0004 


' B 

4(15) 


00G-4R 

El 20 
0 0 7 4 R 

NO GO 

SVC 

: 2 , UNPAK 

UNPACK &. 

00G8R 

El 20 
O078R 


SVC 

2 , L l S i 

PRINT STAT 

00GCR 

El 20 


SVC. 

2 , PAUSE 

UA.IT 


0094 R 





0070R 

4300 
0 0 3 0 R 


B 

GO 

TRY AGAIN 

■ 0074R 

0006 

UNPAK 

DC 

6, MESS • 


/; ■- ■ r. 

00S4R 





007 SR 

0007 

0000 

492-F 

LIST • 

DC 

7 4 12 ,0‘T/0 ERR 



4F20 

4552 

5220 




* 

008 4 R 


NESS 

DS 

i 4 


008 SR 


■ REG 

DS 

..12 


009 4 R 

0001 

PAUSE 

• DC 

■ 1 


0096R 

4 801 

INPUT 

■ ■ DC 

X y 4S01 ',0,0,0 

'LU 1 RB IN 


"■ ' 0000 
0000 
0 000 

009ER END 



page is 


PAGE 2 



PAGE 


READ EC KERMAN TAPE INTO' i;FORTRAH 5/74 
NO ERRORS 


* .0 

0 0 1 £ R 

* ECKERM 

0 0 O 0 R 

GO 

003OR 

INPUT 

009SR 

LIST 

0 0 7 8 R 

LOOP 

0O3CR 

MESS 

0084R 

NOGO 

0064R 

OK 

0 0 1 8 R 

PAUSE 

0 O 9 4 R 

REG 

008SR 

UNPAK 

0074R 
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Fast Fourier Transform 


SUBROUTINE FOUR 1 T ( DATA , NN , S, IS ION) 
v 0 DATA 'COMPLEX ARRAV 

C- NN=NUMB COMPLEX VALS 

C S-SIN TBL , DIM=DIM ( DATA ) /2 + 1 

C ISIGN-SIGN IN FOURIER 

0 IF ABSCISIGN) . NE , 1 , TABLES SETUP 

C I F E Q . 0,TABLE5 ONLY 

0 NO NN NORMALIZATION 

D I MENS I ON DATA ( 1 ) , S (1 ) 

N “ NN+NN 

. NN£=NN/£ 

H N £ 1 = N N 2 + 1 
NN22-NN2+2 
NNP2«=NN + 2 
XN = NN 

IF C IABS< ISIGM) . EQ . I ) GO TO 1 
DO 4 1 = 1, NNSi 

... THETA = 3 , 1415926535 'K FLOAT ( 1-1 ) /XN 
A * S< I > -SINCTHETA) 

IF( I SIGN . EQ . 0) RETURN : 

1- ■ ; J = 1 

DO S I = 1 , N ; £ 

IF(I.GEJ) GO TO 2 
TEMPR-DATAC J) 

TEMPI = DATA C J+l ) 

DATA< J) -DATACI) 

DATAC J+l ■) 'DATA C 1 + 1 ) 

DAT A ( I ) = TEMPR 
DATACI+1 ) =TEMPI 
£ M=N/£ 

3- IF(J.LE.M) GO TO 5 

-J = J-fi 
M = M/£ 

IFCM.GE.2) GO TO 3 

5 T=J+M 

MM AX * 2 

G I STEP = MMAX + MMAX . 

- NMAX = N/MMAX 
DO 8 M= 1 , UMAX, £ 

IND = C ( M — 1 ) &NMAX ) /£+l 
I F C I HD . GT . NN2 1 ) GO TO 7 
UI =S( IND ) 

IND = NN22-I ND 
UR=S ( IND) 

GO TO y 

? IND « NNP2-IND 

WI-S(IND) 

I ND = NN22-I ND 
UR--SCIND) 

9 CONTINUE 

IF ( ISIGM , LE . 0) UI=-UI 
DO 8 I»M,N, ISTEP 



' -J = I + MMAX 

TEMPR“UR#DATA( J)-UI*DAT 

TEMPI =UR*DATA ( J+l )+UI #D 
DATAC J) * DATAC D-TEMPR 
DATAC J+l ’■ = D A T A ( I + 1 ) - TEMPI 
DATA ( I ) - DAT A ( I ) +TEMPR 
DATAC I + n -DATA C I + i ) +TEMPI 
MM AX = I STEP 

IF CMMAX . LT . N ) GO TO 6 

RETURN 

END 


FOUR IT 

O024R 

FUNC/'SUB 

FOUR1T 

0574R ' 

FLING VAR 

; q 

0 0 0 0 R 

EXT FUNG 

.p . = ; 

000OR 

EXT FUNG 

DATA 

0 0 2 A R 

FORM PAR 

NN 

002CR 

FORM PAR 

s 

0 0 S E R 

FORM PAR 

ISIGN . 

0030R 

FORM PAR 

N 

0S7CR 

I NT 4 VAR 

NN2 

0 5 3 0 R 

I NT 4 VAR 

NN21 

@ 5 3 8 R 

I NT 4 VAR 

NN22 

©590R 

I NT 4 VAR 

NNP2 

059 JR 

I NT4 VAR 

xn 

059SR 

REAL VAR 

. U ■ 

0000R 

EXT FUNC 

I APS 

0 0 0 0 R 

EXT FUNC 

i; 

012ER 

LABEL •’ 

4. 

00E4R 

LABEL 

i •. ' ■ c 

059CR 

I NT4 VAR 

THETA 

05A0R 

REAL VAR 

FLOAT 

0 0 0 0 R 

EXT FUNC 

SIN 

0000R 

EXT FUNC 

j : 

05B4R 

I NT4 VAR 

5" 

027CR 

LABEL 

2 . ’ ' •• 

022SR 

LABEL 

TEMPR 

05B3R 

REAL VAR 

TEMPI 

0 5 B C R 

REAL VAR 

M 

05C0R 

I NT 4 VAR 

3 

023 AR 

LABEL 

MM AX 

0SC4R 

INT4 VAR 

6 

0 2 h d R 

LABEL 

ISTEP 

05C8R 

INT4 VAR 

NM AX 

05OCR 

I N T A VAR 

2 

0 4EAR 

LABEL 

IND 

05DOR 

I NT4 VAR 

7 

033CR 

LABEL 

UI 

05D-4R 

REAL VAR 

UR 

05D3R 

REAL VAR 

9 

03SCR 

LABEL 


0003 ERRORS 


<E<E 



34 


0001 

1.0002 

0003 

0004 

0005 
000G 
0007 

0003 

0009 

0010 
001 1 
0012 
0013 
00 1 *4 
0015 

0 0 1 b 

0017 

0013 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 
002? 
0023 
0029. 

0030 

0031 

0032 

0033 

0034 

0035 

0036 
003? 
0033 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

004 7 
0043 
00 4 9 

0050 

0051 

0052 


APPENDIX C 

Main Analysis Program 


0004 R 
OQO 4 R 
0 0 0 4 R 
0 0 0 4 R 
0 0 0 4 R 
0004R 
0 0 © 4 R 
© © 1 4 R 
0 0 3 0 R 


0©3ER 
0 0 5 6 R 
006AR 
003 AR 
O 0 9 C R 


•0 THIS PROGRAM IS TO ANALYZE AND DISPLAY EOPAP SEASAT 
0 6H0RT PULSE RADAR DATA. DEPENDING. ON THE FROGRAM 
€ MODE SELECTED ;• PULSES , PULSE SPECTRA OR PEAK 
0 HISTOGRAMS CAN BE DISPLAYED. EITHER THE AVERAGE 
C SPECTRA OR THE AVERAGE OF THE SPECTRA CAN BE 
€ SELECTED. MEANS AND RMS VALUES- CAN BE EVALUATED. 

•C DATA ARE ADJUSTED FOR POWER-. 

C 

IMPLICIT INTEGERS (I-N) 

DIMENSION IPRC168) ^ 

D I MENS ION INP ( 5344 ) , TBL ( 1 024 ) , DAT ( 5© ) , PULbE^dbb ) 
DIME N S 1 0 N XLAB(b) , PLOT (101 ) , H A V ( 2 to 6 ) , ui I N D ( 2 b b .* 
COMPLEX Z ( 1024 ) JCMPLX . 

EQUIVALENCE < Z C 1 ) , I NP ( 1) ) . 

DATA PI , N , N TUP , LO , IHI/3 , 141592? , 128 ; 1024 , 23 , ib©x 
DATA DOT .STAR, BLANK .SLASH, E S , A A , P E E x 1 H . , 1 H %■ , 1 H , 
CALL FOUR 1 T (Z, NTUP , T BL , 0 ) 

C READ I N M 0 DE DESIRED 
101 WRITE C 0,100) 

100 FORMAT C11H ENTER- MODE X i 

READ (0,2) MODE .-• 

H 2 = M X 2 
•N21-N2+1 


1HX, 1HS, iHA, 1 


0 0 A 3 R 



NTUP21 “NTUP/2+1 : . 


0©BER 



Nl-N+1 


00CAR 



XN = H 


00D8R 



XM- P I XX N 



C 

CALC 

UL ATE WINDOW WEIGHTS 


00E4R 



DO 42 1 = 1, N 


00ECR 

42 


W I N D ( I ) = 1 . — . 8'3 6 3 b 3 6 4 % C 0 S ( F L 0 A T 

( I - 1 ) * X M ) 


C 

GET 

LABEL FOR -SPECTRA 


013ER 



DO 23 K= 1,6 


0 1 4 6 R 

23 


XL AB C K ) » FLOAT (K-l ) 



c 

NOW 

GET INPUT ' PARAMETERS : 


017CR 



WRITE ( 0,4) 


0 1 9 4 R 

4 


FORMAT (6H LABEL ) : ' 


01A2R 



READ (0,5) DAT 


0 i 0 4 R 

S 


FORMAT C 5©Ai ) 


01CER 



UR I TEC ©, 30) 


©1E6R 

. 30 


FORMAT ( 1 3H ANGLE , M I CSEC ) 


01 FOR 



RE AD ( 0 , 31.) ANGLE, DUR )• 


0224R 

31 


FORMAT ( 8F5 . 1 > 


0 2 3 0 R 



WRITE (3, S) DAT 


0252R 

to 


F OR M A T ( I H 1 , 5 0 A !)•• 


0 2 6 0 R 



ANG = ANGLE^P I X 1 80 . 


0270R 



DELTA = 160 . /< DUR5KFL0AT ( NTUP ) ) : 


0288R 



UR I TE ( 0 , 1 ) 


0 2 A 0 R 

. 1 


FORMAT (7H FILE*?,, 7H I THRESH ) 


02BAR 



READ (0.2) IFILE; ITHRSH 


02E2R 

2 


FORMAT (SIS) 


02ECR 



WRITE (3,3) IFILE , ANGLE , DUR 


031 OR 

3 


FORMAT ( 6H FILE , 1 3 , 5X , FS . 1 , 8H 

DEGREES , 5 X , F 6 


MICSEC) 
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0053 

0352R 


. .WRITE (3, 99) ITHRSH 

0054 

0372R 

99 

FORMAT ( 15H THRESHOLD » ,15)5 

0055 

' 038ER 


UR I TEC 0,7) 

0056. 

03A6R 

7 

FORMAT C1SH BLOCKS, SIZE) 

0057 

03BCR 


READC0, 2) IBLKS, ISIZE 

C05S 

0 1 4 R 


' ■ SIZE- ISIZE 

O059. 

03F2R 

■8 

FORMAT C1H , I5,32H AVG SQUARED PULSES , LOG SPECTRUM) 

0O60 

0 4 2 2 R 

41 

•• FORMAT (31H SPECTRA OF AVG=S,AVG SPECTRA-A <) 

0061 


€ 

BLOCKS ARE INPUT IN 32 TRACE BLOCKS . 

0062 

044 AR 


KO-UNT-33 

0O63 

0452R 


DO 9 IIII-1-, IBLKS 

0064 

.. 0 4 5 A R 


. ■ WRITE <3 ,10)1111- 

0065 

047AR 

10 

' FORMAT (7H BLOCK .,14)' 

0066 

04 SCR 


DO- 15 K = LO,IHI 

0067 

- 04S4R 


IL-K+l-LO 

0068- 

04A4R 


PAVCID-0. 

0069 

‘ 04B8R 


IPRCK) =0 

0070 

. 04CAR 

15 

PULSE < K) -0 . 

0071 


C 

WRITE' (3, 22) XL A B 

0072 

04F0R 


CALL E’KSC 1 > . 

0073 

' 04FSR 


C ALL E C K E R M C I M P ) : 

0074 

■ 0 5 A 0 R 


DO 12 1=1, ISIZE 
•• • ' KOUNT=KOUNT+l 

0075 

0508R 


0.076 

O 5 1 4 R 


IF (KGUNT . LE . 32 ) GO TO 13 

007-7 

052SR 


M IN -2TT-1 4 ... 

0078 


c 

GET NEXT INPUT BLOCK OP 32 TRACES 

0079 . 

. 053AR 


CALL ECKERM C IMP > 

0080 

y 54 2R 


DO 11 I 11=1 , 32 

0081 

- 0 5 4 A R 


I J= < 1 1 1-1 ) *167 

0082 

05SAR 


"DO 11 K = L 0 , I HI 

0083 

© 5 6 2 R 


J-IJ + K 

0084 

. 05SER 

11 

NIN = MIN0 ( MIN, INPC J) ) 

0085 

05B4R 


:KOUNT=l 

0086 

05BCR 

13 

I J = C K 0 U N T - 1 ) X 1 67 . 

0087 . 

05 OCR 


DO 14 K = LO , I HI 

0088 

05D4R 


•J-IJ+K 

0089 

05E0R 


X = I H P ( J ) - M I N 

0090 

05FCR 


XX -X SIX 

0091 


■ 0 

KL=K+l-LO 

0092 


c 

ZCKL) =CMPLXC XX , 0 . ) TUIND ( KL ) 

0093 

06©8P 

14 

PULSE (K ) -PULSE CK > +XX 

0094 

' 063ER 


IFCMQDE . N£ . 1 )GO TO 12 ' - 

0095 


•C 

CALCULATE PEAK 1 . HISTOGRAM FOR VALUES ABOVE THRESHOLD , 

0096 

‘ 0650R 


DO 70 K=LO,IHI 

0097 

065SR 


J-IJ+K 

009S 

0664R 


IFC INPC J) . LT . ITHRSH) GO TO 70 

0099 

0680R 


I F < INP C J ) . GT , INPC J-l ) . AND. INPC J ) . GT .INPC 1+1 ) ) IPRCK > 

0100 

' O6F0R 

70 

CONTINUE 

0101 

O702R 


GO TO 12 

0102 


C 

GET PULSE FOURIER TRANSFORM 

0103 

0 7 0 6 R 


CALL FOUR IT C Z, N, TBL, 1 ) 

0104 

0 7 1 4 R 


DO 40 K = 1 , N 2 1 


IPR ( K ) + 1 
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0105 

©7 1 OR 

40 

' © 1 0 6 

0 7 6 A R 

12 : 

0107 

077CR 


010S 

0734 R 

62 

© 1 09 

07AAR 


©110 

07BCR 


0111 

07DCR 

74 

0112 


C Of 

0 1 1 3 

07FCR 


0114 

0304R 


0115 

080CR 


0116 

0 3 i 4 R 


0117 

081CR 


0118 

0824R 


0119 

032 OR 


0120 

083 AR 


0121 

© 8 6 0 R 


©122 

0886R 


0123 

0 8 9 E R 


0124 

08 AAR 


0125 

08BAR 


0 i 28 

08CER 

97 

0127 

08E0R 


0128 

© 3 1 . U K 


0129 

0SF8R 


0130 

©9 1 SR 


0131 • 

- 0940R 

98 

0132 

0 9 6 A h! 


0133 

096ER 


0 i 3 4 

“ ©99ER 


0T3S 

09DER 


0136 

09FER 

96 ■ 

0137 

©A1CR 


0138 

■ 0A42R 

95 

0139 

0 A 58 R 

102 

0140 

0A58R 


0141 

0A6OR 


©142 

0 A 9 0 R 


0143 

0AA4R 


014 .4 ' 

I 0AE8R 

75 

0145 

©AFCR 


0146 

©B 10R 

73 

0147 

0B22R 


0148 


C N 

0149 

0B26R 

72 

0150 

0B2ER 


0151 

0B36R 


0152 

0B3ER 


0153 

0B46R 


0154 

0 B b H R 


0155 

0B7ER 

64 

0156 

0BB2R 



PAV ( K > -PAV ( K ) + CABS C 2 C K) ) 

CONTINUE 

■DO' 62 1*1,101 ’ ■ 

PLOT < I ) = BLANK 
IF < MODE . HE , 1) GO TO 72 
WRITE ( 3 , 74 ) J.SI2E 

FORMATC1H ,15, 15'H PEAK HISTOGRAM ) 

/ U L A T E MEAN;, RMS FOR HISTOGRAM 
S UM MR M = 0 . 

ON T R 0 D = 0 . 

SUMSQ= 0 . 

T M I N = 5 © 0 0 . 

T MA X = 0 . 

DO .97 - I =LO , I HI 
XK = I 

IF( I PRC I ) . GT . 0) 

IF ( IPR C I ) . GT , © ) 

X* IPRC I > 

:S U M N R M = S U M N R M + X 
:C NT.R 0 D = C N T R 0 D + XX X R 
S U M F Q « S U M S Q + X X XKX X K 
CONTINUE 

CNTRGD = CHTROD/SU M H R M 
S U ii S 0 = S U M S 0 / SUM N R M 
.R MS = SORT C SU M SO - C H TROD t'C N T R 0 D ) • 

UIR I TE ( 3 , 9S ) 0 NT RO D , RMS 
FORMAT Cl 3H CENTROID * 

GO TO i©2 

■SI-GZ-RMSKRMS- C < TMAX-TMIN-1 . )XXc±;/i.2 
SI G2P * DURX IS© . XSQRT ( ABS ( SIGZO ) 7 ( IS© 
WRITE <3, 96)5 IGZP 

FORMAT ( 1 5H S I GMA 2 - /Fid.bJ 

0. ) WRITE (3, 95) 

NEGATIVE ’ ) 


TM I N = AMT N 1 ( TM I N , XK ) 
TM AX - A MAX 1 ( TMAX , XK ) 


F S' . 0 , 7 H RMS 


IF ( SIGZ . LT 
FORMAT ( 1 2H 
CONTINUE 
DO 73 K-LO.IHI 


+ . 5 


IX- FLOAT C IPRC K 3 ) /SIZE* TO© 

PLOT (IX) -PEE 

WRITE C 3 , 75 ) K , I PR ( K ) , PLOT 
■FORMAT (1H , 215, 101 Al) 

PLOT ( IX ) = BLANK 
CONTINUE 
GO TO 76 

PLOT PULSE' AVG, SPECTRUM 
XMAX = © . 

CNTROD = 0 . 

SUMNRM = 0 . 

DO 64 K = LO , I HI 

CNTROD = CMTROD+PULSE C K ) XFLOAT ( K ) 
SUMNRM = SUMNRM + PULSE ( K ) 
XMAX-AMAX1 C XM-AX , PULSE C K ) ) 

CNTROD = CN T ROD/ SU M NRM 


= , FS . 0 > 

XCOSC ANG) ) 
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37 


0157 

0BBER 


0158 

0BD0R 


0159 

0BF0R 


v 0 1 GO 

OC10R 

65 

0161 

0C2AR 


0 1 fo 2 

0C3SR 


0163 

0C3ER 


0164 

0C4ER . 


0165 

0 C 7 0 R 


0166 

0C34R 


0167 

0 c A £ H 

61 

0163 

0CC0R 


0169 

CCD4R 

16 

0170 

0D1CR 


0171 : 

0 D 2 4 R 

60 

0172 

0D54R 


0173 


C 

017.4 

0 D 6 S R 

200 

0175 

0 D 6 E R 


0176 

0.D90P 

22 

0177 

. 0DBER 


0173 



0179 

0DC6R 


0120 

0 D E 2 R 


01 31 

0 E 0 S P 


0132 

0E12R 


0183 

0E2SR 


0184 


C -I 

0185 


c 

0186 


c 

0.1.87 

0E54R 


0138 

0E60R 


0189 

0E63R 


0190 

0E7CR 


0191 

0EB6R 

19 

0192 . 


c 

0193 

0ECSR 


0194 

0EDCR 


0195 

OF0ER 

2 0 

0136 

OF2ER 

13 

01-97 

0F40R 

76 ■ 

0193 

0F58R 

21 

01-99 

0F62R 

9 

0200 


C NO 

020 1 . 

0F74R 


0202 

0F7SR 



15) 


. U 
IPR 
I HP 
Z . 
TBL 
DAT 


0000R 

0F7CR 

10CCP 

10CCR 

3A8CR 

“4A8CR 


EXT 

INT2 

I NTS 

CMPX 

REAL 

REAL 


I C H TR = C H T R 0 D + , 5 
UR ITE (3,8) ISIZ.E 
U R I TE C 3 , 6 5 ) I C N T R 
F OR MAT ( i 3H CENTROID 
XNORM= 100 . / X M A X <: 

DO .16 K - 1 , N 
I J - LO + K- 1 

IX -PULSE ( I J ) % X N 0:R M + , 5 
PLOT C IX ) = PEE 
WRITE (3,61)10, PLOT 
F O R M A T C 1 H 13 1 @ i A 1 ) 

'PLOT ( IX) - BLANK 

2 ( K ) = C M P L X C P U L S E ( I J ) /S 1 2 E , 0 . ) 

DO 60 K = N 1 , MTUP 
2 C K. ) = C MP L X C @ . , 0 , ) 

CALL^FOURIT ( 2,- NTUP , TBL, 1 ) 

U R I TE ( 2 , 200) < 2 ( K K ) , K K = 1 , N ) 

FORMAT < 4ES0 , 9 ) 

U R I T E ( 3 , S 2 ) XL A B 
FOR MA T ( 1 0 H F R E QU'ENO y 
DO la K=1 , NTUP2 I- 
P A V C K. ) = P A V C K ) /SIZE 
D = F L 0 A T ( K-l ) 5KDE.LT A 
X= ; CABS ( 2 ( K ) ) y'X N . 

X=ALOG10(X) 

I X = 2 O . X X + . 5 
I> : >MIN0 ( MAX0 ( 1 , IX + i ) 

V * AL.0 G 1 0 C P A V ( K ) ) 

rv = 20 , XY + S I GN < . 5 , Y) 

IV* M I N 0 ( ri A X 0 ( 1 
PLOT ( 1 ) -SLASH 
DO 19 1=2 101 
PLOT ( I ) * BLANK 

I F C ( C I - i ) / 4 ) X 4 . EQ , I - 1 ) PLOT ( I ) = DO T 

CONTINUE 

PLOT ( I Y ) = AA 

PLOT C IX ) =ES 

■WRITE* 3, 20) D , X , PLOT 

FORMAT ( 1H , FS . 2,4X, F7 . 4 SX 101A1) 

CONTINUE ■ •’ 

WRITE* 3 , 21 ) 

.F OR M A T ( 1 H 1 ) 

CONTINUE 
GET NEXT INPMT 
GO TO 101 
END : 

FUNG 

V A R : 

VAR 
VAR 
VAR 
VAR 


4;< IKS , 3 X „ 5* F4 . 1 1 SX 5 , F4 . 1 ) 


101 ) 


IY+1 ) , 101 ) 


f SSSSSS 
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' PULSE 4B54R 
X' PE 4FS4R 
PLOT 4F6CR 
pAV 5100 R 
WIND 5500R 
CMPLX 00O0R 
p-I ■ 00OSR 
N 000CR 

■ NTUP 000ER 
LQ 001 0 R 

I HI 0012R 
DOT 0014R 
-STAR 001 SR 
'•BLANK 00 ICR 
- - SLASH 0020R. 
PS 0024R 

A A 002 SR 

PEE 002CR 
1. FOUR IT 0000R 


■ 101 

003ER 

• ' 100 

0 0 B S R 

<5H 

0000R 

2 ! ' ' 

02E2R 

: • MODE 

5 9 0 4 R 

: N2 

590GR 

. N21 

590CR 

■ . NTUP21 

5912R 

HI 

5914R 

XH , ■ 

5916R 

•• . U 

0000R 

XM ? 

591 AR 

42 

00ECR 

i: 

591ER 

. . cos 

0000R 

FLOAT 

0000R 

. 23 . 

014SR 

• K~ -r 

593SR 

• 4 .V ■’ 

0194R 

&: . 

01C4R 

30 

01E6R 

31 

0224R 

ANGLE 

593 AR 

■ DUR 

593ER 

6 

0252R 

ANG 

5942R 

- DELTA 

594 AR 

!■' 

02A0R 

- IFILE 

5952R 

-- ITHRSH 

5954R 

3 

031CR 

99 

0372R 

7. 

0 3 A 6 R 


REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
EXT FUNC 
REAL VAR 
I NTS VAR 
I NTS VAR 
INT2 VAR 
I NT 2 VAR 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
EXT FUHG 
LABEL 
LABEL 
EXT FUNC 
LABEL 
I NT 2 VAR 
INT2 VAR 
I NT 2 VAR 
I NT 2 VAR 
INT2 VAR 
REAL VAR 
EXT FUNC 
REAL VAR 
LABEL 
INT2 VAR 
EXT FUNC 
EXT FUNC 
LABEL 
I NT 2 VAR 
LABEL 
LABEL 
LABEL 
LABEL 
REAL VAR 
REAL VAR 
LABEL 
REAL VAR 
REAL VAR 
LABEL 
INTE VAR 
I N T 2 VAR 
LABEL 
LABEL 
LABEL 



iblks 
ISI2E 
' SIZE 
8 

41' 

• KOUNT 

9 

I II I 

10 
IS 
IL 
BKS 

eckerm 

IS 

13 
MIH 
.1 
11 

• Ill 
IJ 

' J, 

- H I N 0 

14 

X- ■- 

XX 

70 

40 
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s a 
75 
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0 3 F S R 

04S2R 
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0FSSR 

596 4 R 

0 4 7 A R 
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0O00R 

076 AR 

05 BCR 

5970R 

000AR 

056ER 

S376R 

S97SP 

597ER 

0000R 

0603 R 

5980R 

5 9 S 4 R 

06F 0R 

07 1 CP 

0000R 

0784R 

0B2SR 


74 07DCR 

' SUMNRM 598SR 


■ CNTROD 

5 93 OR 

■ SUMSQ 

5990R 

TMIN 

5994P 

TMAX 

599CP 

97 

0SCER 

XK 

59A0R 

AMIN1 

0O00R 

AN AX 1 

0000R 

RNS 

5 9 A 4 R 

SORT 
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93 

0S40R 

102 

0A53R 

SIGZ 

S5A8R 

. R 

0000R 

SIGZP 

59B0R 

ABS 

0000R 

96 

09FER 

95 

0A4SR 

73 

0B10R 

IX 

59C0R 


I NTS VAR 
INTS VAR 
REAL VAR 
LABEL 
LABEL 
INTS VAR 
LABEL 
INTS VAR 
LABEL 
LABEL 
INTS VAR 
EXT FUNC 
EXT FUNC 
LABEL 
LABEL • 
INTS VAR 
EXT FIJNC 
LABEL 
INTS VAR 
INTS VAR 
INTS VAR 
EXT FUNC 
LABEL 
REAL VAR 
REAL VAR 
LABEL 
LABEL 
EXT FUNC 
LABEL 
LABEL • 
LABEL • 
REAL .VAR 
RtAL VAR 
REAL VAR 
REAL VAR 
REAL VAR 
LABEL 
REAL VAR 
EXT FUNC 
EXT FUNC 
REAL VAR 
EXT FUNC 
LABEL 
LABEL 
REAL- VAR 
EXT FUNC 
REAL VAR 
EXT FUNC 
LABEL 
LABEL 
LABEL 
INTS VAR 
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. V 

75 

76 

XHAX 

6-1 
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£5 

XNORM 
1 6 
G 1 
SP 
GO 
20 O 
22 
18 
D 

ALOG10 
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19 

20 
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0D62R 
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0F2ER 

B9D4R 

0000R 
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0000R 


EXT FUNC 
LABEL 

label 

REAL VAR 


label 
1 N T 2 VAR 
label 

REAL VAR 
LAB.EL 
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EXT FUNC 
LABEL 

label 

label 

label 

PEAL VAR 
ext FUNC 
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LABEL 
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label 
ext FUNC 


0000 ERRORS 
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